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Abstract 

We consider the Langevin equation with multiplicative noise term which 
depends on time and space. The corresponding Fokker-Planck equation in 
Stratonovich approach is investigated. Its formal solution is obtained for 
an arbitrary multiplicative noise term given by g(x,t) = D(x)T(t), and the 
behaviors of probability distributions, for some specific functions of D(x), are 
analyzed. In particular, for D{x) ~ l^l -6 ^ 2 , the physical solutions for the 
probability distribution in the Ito, Stratonovich and postpoint discretization 
approaches can be obtained and analyzed. 
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Typeset using REVTgX 



In these last decades, anomalous diffusion properties have been extensively investigated 
by several approaches in order to model different kinds of probability distributions such as 
long-range spatial or temporal correlations [I]. For instance, the well-known cases are the 
Langevin and the corresponding Fokker-Planck equation, and the master equation. The 
other ones we could mention are the generalized Langevin equations [2], the generalized 
Fokker-Planck equation with memory effect [3], generalized thermostatistics [4], generalized 
master equations [5], continuous time random walk [6], and fractional equations [I]. These 
approaches have been used to describe numerous systems in several contexts such as physics, 
hydrology, chemistry and biology. 

The well-established property of the normal diffusion described by the Gaussian dis- 
tribution can be obtained by the usual Fokker-Planck equation with a constant diffusion 
coefficient (without the drift term). Anomalous diffusion regimes can also be obtained by 
the usual Fokker-Planck equation, however, they arise from variable diffusion coefficient 
which depends on time and/or space. On the other hand, in the view of Langevin approach 
it is associated with a multiplicative noise term. In other approaches such as the generalized 
Fokker-Planck equation (nonlinear) and fractional equations, they can describe anomalous 
diffusion regimes with a constant diffusion coefficient. 

In this work, we investigate the Fokker-Planck equation with variable diffusion coefficient 
in time and space, in the Stratonovich approach. We show that for a multiplicative noise 
term separable in time and space, g(x, t) = D(x)T(t), we can obtain a formal solution for the 
probability distribution. We also analyze the behaviors of probability distributions for some 
specific functions of D(x) which can manifest interesting properties such as non-Gaussian 
distribution, combination of behaviors like Gaussian (for small distance) and exponential 
(for large distance), and combination of behaviors like Gaussian (for small distance) and 
power law decay for long distance. Also, we can obtain many bimodal distributions for 
different forms of D(x). 

Now, we consider the following Langevin equation 



where & is a stochastic variable and Tj(t) is the Langevin force. We assume that the averages 
(Vi(t)) = and (r^r^f)) = 2<J <i <5(t — *) [3]. For h = and g = D, Eq.(l) describes the 
Wiener process and the corresponding probability distribution is described by a Gaussian 
function. In the case of <?(£, t) variable, some specific functions have been employed to study, 
for instance, turbulent flows (g(x,t) ~ |x| a t 6 ) [7], relativistic Brownian motion (h(p,t) ~ p 
and g(p,t) ~ (1 + (p/mc) 2 ) 1 ^) [8] and classical harmonic oscillator [9]. By applying the 
Stratonovich approach in a one-dimensional space [3], we obtain the following dynamical 
equation for the probability distribution 



(i) 



dW(x,t) 
dt 



d[h(x,t)W(x,t)\ d 



dx dx 



d(g(x,t)W(x,t)) 



dx 



(2) 



Hereafter we consider the drift term equal to zero (h — 0). We also consider that the 
multiplicative noise term g(x,t) is separable in time and space, g(x,t) = D(x)T(x). Then, 
we can rewrite Eq. (2) as the following manner 



with 



p(x,t) = D(x)W(x,t) , 



(4) 



and 



di 
di 



dx 



T\t) 



1 



dx D(x] 
Eq. (3) has the following formal solution, 

exp 

p(x, t) = A 



(5) 



(6) 



it 



(7) 



where A is a normalization factor. We note that for D{x) = \f~D and T(t) = 1 we recover 
the Wiener process. 

We can now explore some spatial features of the probability distribution of solution (4), 
for some specific forms of D(x): 
First case - We consider 



x 



x 



VD{1 + a\x\ c ) b ' 

where a is a positive real number. From Eqs. (4), (6) and (7) we obtain 

[1 + a(l — be) \x\ c ] exp 



(8) 



W{x,t) = Ci- 



ADi(l+a\x\ c ) 2b 



Vra + a\x\ c ) 1+b 



(9) 



where be < 1 in order to maintain W(x,t) positive. For a = or b = we recover the 
Wiener process. 

In particular, for c = 1/6, we have 



exp 



W{x,t) = C 2 - 



4Dt(l+a\x\ 1/b ) 



2F 



Vi(l + a 



x 



i/by+ b • 



(10) 



In this process, the behavior, for small a and x, is like a Gaussian function, whereas for large 
distance, the exponential term converges to a constant value. Therefore, for large distance, 

/ \l/b\ 

the dominant term is the multiplicative factor 1/ ( 1 + a \x\ 1 which approximates to the 

asymptotic power law x^ 1+b ^ b . We note that, for instance, the power lay decay is present in 
the fractional and nonlinear approaches [1,10]. 
For 6 = 1 and c = 1/2 we obtain 
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W(x,t) = C 3 - 



2 + a^j\x\ exp 


X 2 


ADt 


^l+a^/\x~^ 


■1 




x\) 



(11) 



In this process, the behavior, for small a and x, is like a Gaussian function. For large 
distance, we have a exponential decay basically. We note that the exponential decay has 
been observed in pair dispersion in two-dimensional turbulence [11]. 

Moreover, for be < 1/2, the decay of the solution (9) is essentially compressed Gaussian 
shape, whereas for 1/2 < be < 1, the decay is essentially stretched Gaussian shape. It is 
interesting to emphasize that the solution (9) can have a similar asymptotic non-Gaussian 
shape of the random walk model and time-fractional dynamic equation [1]. The asymptotic 
shape of the random walk model and time-fractional dynamic equation is given by 



W(x,t) ~ C 4 r~C~ exp -C 5 £ 



(12) 



where £ = \x\ jt a l 2 . This shape can be obtained from the solution (9), for large distance, by 
taking be = (1 - a)/ (2 -a),t = t a /( 2 " a ) and T 2 (t) = _a_f2(a-i)/(2-a)_ 
Second case - We consider 



D(x) = VD \x\ 



(13) 



where 9 is a real parameter. We should note that the diffusion coefficient (13) has been used 
to describe the diffusive process on a fractal [12]. The probability distribution (4) for the 
spatial multiplicative noise term (13) is given by 



W(x,t) 



e 


\x\ 2+e 


\x\ 2 exp 


D(2+e) 2 t 



(14) 



In this process, we have the bimodal states. In fact, we can construct many bimodal states 
by choosing different functions for D(x). For 9 = we recover the Wiener process. Basically, 
for large distance, the probability distribution (14) has a non-Gaussian decay. The second 
moment related to this process is given by 



x 



[D 2 (2 + 9Y]~ e T 



6+0 



2(2+6>) 



2 

t 2+8 



(15) 



2/(2 -a), 
~ t a which 



The solution (14) also reproduces the asymptotic shape (12) by taking 2 + = 
t = r/( 2 - Q ) and T 2 {t) = ^t 2 ( a - 1 )/( 2 ~ a ). The second moment (15) yields (x 2 ) 
corresponds to the same behavior of the time-fractional diffusion equation [1]. For this 
process, the multiplicative noise term corresponds to (g(x,t) ~ |x| a t b ) which has the same 
form suggested by Hentshel and Procaccia to study the turbulent system [7]. 

We can now compare with the solution obtained by the Ito approach for the same 
Langevin equation using (13). The solution has been obtained in [13], and for T(t) = 1 
it is given by 
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The second moment yields 

„ t^e . (17) 

We see that these two approaches give different behaviors for the probability distribution 
due to the multiplicative factors. However, the second moment of these two approaches give 
the same behavior. Further, both the distributions present the bimodal states. It is also 
interesting to compare with an other approach which uses the postpoint discretization rule 
[8,13]. For this last case, the probability distribution does not present the bimodal states, 
however, its second moment has the same power law behavior of the Ito and Stratonovich 
approaches. We see that three different approaches give different behaviors, but they give 
the same power law behavior for the second moment. 

In summary, we have investigated the usual one-dimensional Fokker-Planck equation 
with variable diffusion coefficient in the Stratonovich approach. We have considered a very 
general class of the multiplicative noise term, g(x, t) = D(x)T(t), and we have presented the 
formal solution for the probability distribution. Using the formal solution, we have analyzed 
some particular solutions by choosing simple functions for D(x). We have shown interest- 
ing behaviors for the probability distribution such as non-Gaussian, exponential and power 
law decays for large distance. As we can note that the usual Fokker-Planck equation can 
describe many different anomalous processes with many different behaviors. The introduc- 
tion of a time dependent multiplicative noise term may be necessary for the cases of more 
complex systems such as turbulent systems, as suggested by several authors [7]. In fact, 
we have shown that the asymptotic shape of the random walk model and time-fractional 
dynamic equation can be obtained from the solutions described in this work with the time 
dependent multiplicative noise term. Further, we have also shown that the solutions of 
the Ito, Stratonovich and postpoint discretization approaches, for D(x) = \f~D \x\ ~ e l 2 ^ de- 
scribe different behaviors, but their second moments describe the same behavior. If the 
diffusion coefficient D(x) = \f~D \x\ 9 ^ 2 may describe exactly a real physical system by the 
Fokker-Planck equation, then further information of the microscopic structure of the system 
is necessary in order to choose which of the above approaches is the correct one or, simply, 
which of the above approaches can fit the experimental data. 
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